AD-A244  098 


MEMORANDUM  REPORT  BRL-MR-3953 


A  COMPUTER  PROGRAM  TO  REDUCE 
DIGITIZED  YAW-CARD  DATA  FROM 
FINNED  KINETIC-ENERGY  PROJECTILES 


RICHARD  A.  PENNEKAMP 


DECEMBER  1991 


APPROVED  FOR  PUBUC  RELEASE;  DISTRIBUTION  IS  UNLIMITED. 


U.S.  ARMY  LABORATORY  COMMAND 


BALLISTIC  RESEARCH  LABORATORY 
ABERDEEN  PROVING  GROUND,  MARYLAND 


electk 

■JAN  a  8 1992 


NOTICES 


Destroy  this  report  when  it  is  no  longer  needed.  DO  NOT  return  it  to  the  originator. 


Additional  copies  of  this  report  may  be  obtained  from  the  National  Technical  Information  Service, 
U.S.  Department  of  Commerce,  5285  Port  Royal  Road,  Springfield,  VA  22161. 


The  findings  of  this  report  are  not  to  be  construed  as  an  official  Department  of  the  Army  position, 
unless  so  designated  by  other  authorized  documents. 


The  use  of  trade  names  or  manufacturers'  names  in  this  report  does  not  constitute  indorsement 
of  any  commercial  product. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
0MB  No.  0704-0188 


Public  feDOrting  burden  for  this  collection  of  information  is  estimated  to  average  i  nour  oer  'ew>orse.  including  the  time  tor  reviewing  instructions,  searcning  enstmg  data  sources, 
gathering  and  maintaining  the  data  needed,  and  como^eting  and  reviewing  the  collection  of  information  Send  comments  regarding  this  burden  estimate  or  any  other  asoect  of  this 
colleaion  of  information,  including  suggestions  for  reducing  this  ouroen.  to  ^Washington  neadouaners  Services.  Directorate  for  information  Ooerations  and  Reports.  1215  Jefferson 
Oavis  Highway.  Suite  t204,  Arlington.  VA  22202<4302.  and  to  the  Office  of  Management  and  Budget  PaperworK  Reduaion  Protect  (0704-0188).  Washington.  DC  20S03. 


1.  AGENCY  USE  ONLY  {Leave  blank) 


2.  REPORT  DATE 

Deceaber  1991 


3.  REPORT  TYPE  ANO  OATES  COVERED 

Final.  Apr-May  91. 


4.  TITLE  AND  SUBTITLE 

S.  FUNDING  NUMBERS 

A  Computer  Program  to  Reduce  Digitized  Yaw-Card  Data 
from  Finned  Rlnetlc-Energy  Projectiles 

6.  AUTHOR(S) 

1L162618AH80 

Richard  A.  Pennekamp 

7.  PERFORMING  ORGANIZATION  NAME(S)  ANO  AOORESS(ES) 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

9.  SPONSORING /MONITORING  AGENCY  NAME(S)  ANO  AODRESS{ES) 

Ballistic  Research  Laboratory 

10.  SPONSORING /MONITORING 

AGENCY  REPORT  NUMBER 

ATTN:  SLCBR-DD-T 

Aberdeen  Proving  Ground,  MD  21005-5066 

SEL-MR-3953 

11.  SUPPLEMENTARY  NOTES 

12a.  DISTRIBUTION /AVAILABILITY  STATEMENT 

12b.  DISTRIBUTION  CODE 

Approved  for  public  release;  distribution  Is  unlimited. 

13.  ABSTRACT  {Maximum  200  worth) 

A  method  of  utilizing  an  electronic  digitizer  to  reduce  yav  card  data  from 
finned,  kinetic  projectiles  has  been  developed.  A  listing  of  a  FORTRAN  program 
that  uses  this  digitized  data  Is  provided. 


14.  SUBJECT  TERMS 


Yaw 

Computer  Programs 
Kinetic  Energy  Projectiles 


17.  SECURITY  CLASSIFICATION 
OF  REPORT 

UNCLASSIFIED 


NSN  7540-01-280-5500 


18.  SECURITY  CLASSIFICATION 
OF  THIS  PAGE 

UNCLASSIFIED 


19.  SECURITY  CLASSIFICATION 
OF  ABSTRACT 


15.  NUMBER  OF  PAGES 

30 


16.  PRICE  CODE 


20.  LIMITATION  OF  ABSTRACT 


UNCLASSIFIED 


UNCLASSIFIED 


Standard  Form  298  (Rev  2-89) 

Pre»crfbed  by  ANSI  Sfd  /39-«8 
298-102 


Intentionally  left  blank. 


ii 


TABLE  OF  CONTENTS 


LIST  OF  FIGURES .  v 

1.  INTRODUCTION .  1 

2.  THEORY .  2 

3.  CONCLUSION .  6 

4.  REFERENCES .  7 

APPENDIX  A:  LISTING  OF  A  FORTRAN  77  PROGRAM  THAT  REDUCES 
DIGITIZED  YAW-CARD  DATA .  9 

APPENDIX  B:  INSTRUCTIONS  FOR  USING  THE  YAW-CARD  REDUC¬ 
TION  PROGRAM .  17 

APPENDIX  C:  HOW  THE  4  LENGTH  IS  DETERMINED .  25 

LIST  OF  SYMBOLS  AND  ABBREVIATIONS .  31 

DISTRIBUTION  LIST  .  33 


iii 


ioosssion  For 

NTIS  GRAM 

DTIC  TAB 

□ 

Unf.nnoiiiiced 

□ 

Just,  1  r  i  c'it  1  or _ 

! 

By. 

^Distribution/ 

Avnllr.bll Ity  C.- 

des 

Dlst 


ft 


.1 


jAvall  auO/or 
SpeolaX 


Intentionally  left  blank. 


iv 


LIST  OF  FIGURES 

Figure  Page 

1  Yaw- Card  Coordinate  System .  1 

2  Relationship  between  the  Yaw-Card  Hole  and  S .  2 

3  Digitizer  Coordinate  System .  4 

Bl  Projectile  Information  Required  to  Execute  Program .  19 

B2  Digitizing  Pad  Layout  .  20 

B3  Required  Digitized  Data  to  Execute  Program  .  23 

Cl  Case  1:  Fin  Blades  in  Yaw  Plane .  29 

C2  Case  2:  Fin  Blades  Perpendicular  to  Yaw  Plane .  30 


V 


1.  INTRODUCTION 


Witness  panels,  commonly  referred  to  as  yaw  cards,  are  used  in  ballistic  testing  to  de¬ 
termine  the  pitch  (a)  and  yaw  (/?)  of  a  projectile  as  a  function  of  range.  Yaw  cards  are 
traditionally  mecisured  by  a  device  called  a  yaw-card  template.  This  device  is  calibrated  to 
measure  the  total  angle  of  attack,  6,  where  8  is  related  to  a  and  /?  : 


a  =  Ssind 

(1) 

0  =  ScosO 

(2) 

The  angle  0  is  defined  in  Figure  1.  In  order  to  obtain  6,  an  instrument  that  measures  angles 
(usually  a  protractor)  must  be  utilized.  The  a  and  /?  data  are  usually  input  to  some  computer 
program:  for  example,  to  determine  flight  characteristics  of  a  projectile. 


•pS  SEEN  FROM  THE  MUZZLE" 


Figure  1.  Yaw-Card  Coordinate  System. 

The  aforementioned  technique  for  measuring  yaw-cards  is  adequate,  but  requires  manu¬ 
ally  entering  data  into  a  computer  after  measuring  two  quantities  with  two  different  tools. 
A  method  that  uses  an  electronic  digitizer  h«is  been  developed  to  quicken  the  data-reduction 
process.  It  is  a  procedure  that  acquires  data  that  can  be  used  to  calculate  q,  /3,  and  the 
fin  orientation  angle  0.  The  V’  value  is  used  to  determine  the  spin  rate  ais  a  function  of 
range  (Pennekamp  1990).  The  data  obtained  by  the  digitizer  can  either  be  processed  by  the 
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digitizer’s  computer  or  transferred  to  anotlier  computer  for  processing.  Tiiis  eliminates  the 
necessity  of  manually  entering  the  data  into  the  computer. 

This  report  describes  how  the  digitized  data  are  manii)ulaled  to  obtain  the  required 
results.  Appendix  A  has  a  listing  of  a  FORTRAN  program  that  aids  in  the  reduction  of  the 
yaw-card  data.  An  example  of  the  use  of  the  the  program  is  given  in  Appendix  B. 

2.  THEORY 

In  general,  yaw  cards  are  used  to  determine  the  S  value  of  a  projectile  at  a  given  range 
location.  The  8  value  has  the  following  relationship  with  the  hole  in  the  yaw  card: 

Clen  =  ^  (^Os8  -f  ip  sinS,  (3) 

where  C/en  is  defined  in  Figure  1,  r  is  the  radius  of  the  flight  body,  and  is  the  length  of 
the  cylindrical  portion  of  the  flight  projectile  (see  Figure  2).  Aj^pendi.N  C  describes  how  ip 
is  determined.  If  small  angles  (<  10*’)  arc  assumed,  Ec|uation  3  simplifies  to 

f  ..if  C/en  -  r\ 

It  should  be  noted  that  Oc  >  8  for  Equation  4  to  be  valid. 


Figure  2.  Relationship  between  the  Yaw-Card  Hole  and  8. 

The  only  unknown  value  in  Equation  4  is  C/e„.  As  seen  in  Figure  3,  Cun  is  the  distance 
between  PI  and  P2.  If  the  coordinates  of  these  two  points  were  known.  Cun  could  be 
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determined  by 


Cuu  =  ^/{PU-P2,y  +  {P\,-P■2y)\  (5) 

where  the  subscripts  denote  either  the  “x”  or  “y”  value  in  the  traditional  Cartesian  coordinate 
system. 

A  digitizer  is  a  machine  that  provides  the  coordinates  of  a  user-selected  point.  The  user 
places  a  component  of  the  digitizer  on  the  point  and  presses  a  button.  A  signal  is  sent  to  the 
computer  portion  of  the  digitizer  where  the  coordinates  of  the  point  are  determined.  This 
process  will  be  referred  to  as  “digitizing”  a  point. 

With  the  use  of  the  digitizer,  the  locations  of  PI  and  P2  can  be  determined  in  the  XY 
plane.  The  location  of  P2  is  easy  to  acquire;  just  locate  the  extreme  edge  of  the  yaw  card 
hole  and  digitize  the  point.  PI  cannot  be  determined  in  this  straightforward  manner.  PI 
is  the  center  of  a  circle  whose  diameter  is  the  fin  span.  The  extreme  edges  of  the  fin  blade 
imprints  (points  P3-P8)  lie  on  the  circumference  of  this  circle.  (In  reality  it  is  an  ellipse,  but 
for  small  angles  a  circle  is  a  reasonable  approximation.)  The  points  P3  through  P8  can  be 
digitized,  and  they  must  satisfy  the  equation  of  a  circle  in  the  XY  coordinate  system  (Be3'er 
1981), 

fl/in  =  (rir  -  Pn.y  +  (PI,  -  Pn,)^  (6) 

where  is  the  radius  of  the  fin  span,  and  n  has  the  value  of  3  through  8.  The  only 
unknown  values  in  Equation  6  are  the  x  and  y  values  of  PI,  and  these  can  be  determined 
by  applying  nonlinear  least  squares  (Scarborough  1930)  to  Equation  6.  The  nonlinear  least 
squares  technique  is  iterative  and  requires  initial  guesses  to  start  the  process:  reasonable 
estimates  for  the  x  and  y  values  of  PI,  They  can  be  provided  by  digitizing  a  point  that  is 
near  PI,  for  example  P9  in  Figure  3.  After  a  solution  to  the  nonlinear  least-squares  technique 
is  obtained,  the  location  of  P2  is  known.  This  allows  Cien  to  be  determined  (Equation  5) 
and  from  that  information,  S  can  be  obtained  (Equation  4). 

In  order  to  calculate  a  and  the  angle  0  must  be  determined.  The  angle  0  is  given  by 
the  relation 

0  =  <^  +  9O‘',  (7) 

where  (f>  is  the  angle  between  the  reference  line  and  the  line  that  contains  PI  and  P2  (see 
Figure  3).  The  points  PIO  and  Pll  are  on  the  reference  line  and  their  x  and  y  values  can  be 
obtained  with  the  digitizer.  The  x  and  y  values  of  Pi  and  P2  have  already  been  obtained 
in  order  to  calculate  S.  With  the  x  and  y  values  for  points  PI,  P2,  PlO  and  Pll  known,  the 
direction  cosines,  A  and  B,  can  be  calculated  from  the  following  formulas  (Beyer  1981): 

,  F2,  -  Ph 
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Figure  3.  Digitizer  Coordinate  System. 


P9  _  PI 
r>  ^“V  ^ 

1 

“2-1 


/lo  — 


PIU  -  PIO. 


^  j  / 

“11-10 

=  (.1) 

“11-10 

where  c?2-i  is  the  distance  between  P2  and  Pi  and  can  be  calculated  by  using  Equation  5 
(d2_i=C/en)-  The  dii_io  term  is  defined  in  a  similar  fashion.  With  the  direction  cosines 
known,  an  angle  6'  can  be  determined  by  (Beyer  1981): 

(f>  =  cos  ^{AiA2  +  B\B2)-  (12) 

Only  the  principle  value  of  the  arc  cosine  is  computed;  that  is, 

<<j>'  <  180^  (13) 

If  the  angle  (f)  (as  shown  in  Figure  3)  is  less  than  180“’,  then 


and  if  not,  then 


4>  =  <i>'\ 


4>  =  360'’  -  4>' 


A  method  to  determine  if  (f)  is  larger  or  smaller  than  180°  will  now  be  discussed. 

The  first  step  in  determining  if  4>  is  larger  or  smaller  than  180°  requires  that  the  slope  of 
the  reference  line  be  calculated  by  the  following  equation: 

PlQy-PUy 

FlOx-Pllx 

The  slope  from  Equation  16  (Are/)  can  be  infinity  (oo),  zero,  positive  or  negative. 

The  first  case  examined  is  Are/=c>o  (this  occurs  when  P10x  =  Pllr)-  The  p  angle  is  greater 
than  180°  if  Pl^  is  greater  than  P2x.  If  this  is  the  case,  then  Equations  12  and  15  are  used 
to  calculate  (j).  If  not.  Equations  12  and  14  are  used  to  determine  (f>. 

For  all  the  other  possible  values  of  Are/,  the  equation  of  the  line  that  goes  through  the 
Point  PI  and  has  the  slope  of  Are/  must  be  calculated.  This  is  accomplished  by  determining 
the  y-intercept  (6re/)  of  the  line  from  the  following  equation: 

6re/  =  Ply  -  Are/Plx  (17) 


For  Are/  zero  or  positive,  the  line  with  the  slope  of  Are/  and  the  y-iutercept  of  brej  can 
be  used  to  establish  the  following  inequality: 

1/  <  Are/I  +  6re/.  (IS) 
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If  P2’s  X  and  y  values  do  not  satisfy  the  above  inequality,  then  (f>  is  determined  by  using 
Equations  12  and  15.  If  they  do,  then  (j)  is  calculated  by  using  Equations  12  and  14. 

For  Are/  negative,  the  following  inequality  is  utilized: 

?/>  Are/a:  +  6re/.  (19) 

Again,  if  P2’s  x  and  y  values  do  not  satisfy  the  above  inequality,  then  (f>  is  determined  by 
using  Equations  12  and  15.  If  they  do,  then  (f>  is  calculated  by  using  Equations  12  and  14. 

From  this  information,  6  can  be  calculated  from  Equation  7,  and  a  and  ^  can  be  calcu¬ 
lated  by  using  Equations  1  and  2. 

Another  bit  of  information  that  can  be  obtained  from  a  yaw  card  is  the  fin  orientation 
angle,  The  locations  of  all  the  points  (PI,  P3,  PIO  and  Pll)  required  to  determine  il' 
would  have  to  be  obtained  in  order  to  calculate  S  and  0.  With  this  information,  the  direction 
cosines  can  be  calculated,  and  tl’  can  be  determined  in  the  same  fashion  as  (^. 

3.  CONCLUSION 

A  method  of  measuring  the  data  from  a  yaw  card  using  an  electronic  digitizer  has  been 
df  scribed.  P  j'-qui'e:  digitizing  a  rmitc  number  of  points  (~  10)  on  each  card,  and  from 
•hat  i:if  "ma'ion,  tlir  total  angle  of  attack,  the  pitch,  the  yaw,  and  the  fin  orientation  angle 
can  be  obtained.  The  process  produce"-  data  on  computer  files,  which  eliminates  manually 
entering  the  data. 

The  results  from  the  digitized  readings  compare  well  with  the  results  from  reading  tlie 
yaw  cards  with  the  template  and  protractor.  Some  213  yaw  cards  were  read  with  both 
methods,  and  the  average  difference  for  the  S  readings  was  less  tiian  0.05°.  For  the  d’  values, 
the  average  difference  was  0.2°.  (The  0  values  were  not  read  with  the  protractor,  but  their 
results  are  assumed  to  be  similar  to  the  values.)  These  results  suggest  that  no  accuracy  is 
lost  by  using  the  digitizer  to  read  the  yaw  cards. 

Because  there  is  no  loss  ir.  accuracy  and  a  potential  for  time  savings,  the  utilization  of 
a  digitizer  to  aid  in  the  reduction  of  yaw  cards  is  a  logical  stej)  in  tlie  modernization  of 
ballistic-data  reduction. 
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APPENDIX  A: 

LISTING  OF  A  FORTRAN  77  PROGRAM  THAT  REDUCES  DIGITIZED 

YAW-CARD  DATA 
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dimension  x(10),y(10) 
character*72  decrib 

♦File  3  contains  pertinent  projectile  information 

open (3 , f ile= ‘ rdtype ' , status* ' old ' ) 

read(3, 101)decrib 

read(3 ,*)f insp 

read(3,*)plen 

read(3,*)prad 

read(3,*)f inno 

close(3) 

radf in=f insp/2 . 

J  =  1 

******** 

♦File  1  contains  the  digitized  data.  File  2  is  the  output  file 
******** 

open ( 1 , f ile= ’ digdat ’ , status* ’ old ' ) 
open (2 , file* 'output ' , status* 'new' ) 
write (2, 101) decrib 
write(2,103) 

* 

* - 

♦Loop  5  reads  the  data  from  the  DIGDAT  file. 

* 

5  read(l,^,end=6)xr2,yr2 

read(l .♦)xrl ,yrl 
i*l 

2  read(l,^)x(i) ,y(i) 

******** 

♦This  statement  looks  for  the  switch  that  signals  the  end  of  fin- imprints 
♦entries.  It  has  been  arbitrarily  set  at  (3,3). 

4c******* 

if ((x(i) .lt.3) .and. (y(i) .lt.3))gotol 

i»i+l 

goto2 
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numpts=i-l 
read(l,*)xc,yc 
read(l ,*)xdelta,ydelta 
call  circen(x,y ,numpts,xc,yc,radf in) 
pi=2*asin(l . ) 

*phi  and  theta  calculations 

call  angle (xrl.yrl ,xr2,yr2,cal,cbl) 
call  angle (xc,yc,xdelta,ydelta,ca2,cb2) 
phi=acos(cal*ca2+cbl*cb2) 
phi=phi*180/pi 

call  truajig(xc,yc,yr2 ,xr2,yrl ,xrl ,xdelta,ydelta,phi) 
theta=phi+90 . 

if (theta. gt .360 . )theta=theta-360 . 

*psi  calculations 

call  angle(xc,yc,x(l) ,y(l) ,ca2,cb2) 

psi=acos(cal*ca2+cbl*cb2) 

psi=psi*180/pi 

call  truang(xc,yc,yr2,xr2,yrl,xrl,x(l) ,y(l) ,psi) 
psi=psi+90 . 
spaang=360 . /f inno 

4  if ((psi .ge .0 . ) . and. (psi .le . spaang) )goto3 

psi=psi-spaang 
goto4 

♦delta,  alpha  and  beta  calculations 

3  clen=sqrt ( (xc-xdelta)**2+(yc-ydelta)**2) 

delta=asin((clen-prad)/plen) 
delta=delt a* 180/pi 
alpha=delta*sin(theta*pi/180 . ) 
beta=delta*cos (theta*pi/180 . ) 

4>«  4<4I4<4<4<4I 


12 


output  write  statement 


♦ 

write(2,102) j , delta, theta, alpha, beta, psi 
read(l ,*) 
j=j  +  l 
gotoB 

♦ 

* - 

* 

6  close(l) 

close(2) 

101  format (a) 

102  format('  ',i3,5fl2.3) 

103  format('  ' ,//,2x, 'no. ' ,6x, ’delta’ ,5x, 'theta' ,8x, 

:  'alpha' ,7x, 'beta' ,8x, 'psi' ,/) 

end 

★ 

* 

♦This  subroutines  does  the  preliminary  calculations  required  to 
♦determine  the  angle  between  two  lines. 

♦ 

subroutine  angle(xl ,yl ,x2,y2,ca,cb) 

dist=sqrt ((x2-xl)^(x2-xl)+(y2-yl)^(y2-yl)) 

ca=(x2-xl)/dist 

cb=(y2-yl)/dist 

end 

* 

* 

♦This  subroutine  computes  the  center  of  the  fin-blade  imprints. 

* 

subroutine  circen(x,y ,numpts ,x0 ,y0 ,rad) 
double  precision  c(2,2),r(2) 
dimension  x(10),y(10) 
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do  1  i=l,2 
r(i)=0. 
do  1  j=l,2 
1  c(i,j)=0. 

do  4  i=l,numpts 

tl=(2*x0-2*x(i)) 

t2=(2*y0-2*y(i)) 

t3=rad*rad-  ((x(i)-xO)*(x(i)-xO)  +(y(i)-yO)*(y(i)-yO)) 

c(l,l)=c(l,l)+tl*tl 

c(2,2)=c(2.2)+t2*t2 

c(2,l)=c(2,l)+tl*t2 

r(l)=r(l)+t3*tl 

4  r(2)=r(2)+t3*t2 
c(1.2)=c(2.l) 

call  simul(c,r,2) 

count*count+l 

if  (count .gt . 100) return 

if (abs((r(l)) .gt. .001) .or. (abs(r(2)) .gt. .00l))goto5 
return 

5  x0*x0+r(l) 
y0*y0+r(2) 
goto  6 
end 
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* 

^(H^^^^^^iitl**************  **************************************  ************** 
* 

♦This  subroutine  solves  a  system  of  equations. 

* 

SUBROUTINE  SIMUL(C,R,N) 

DOUBLE  PRECISION  TEMP.R(2) ,C(2,2) 

K=1 

9  TEMP=1./C(K,K) 

J=K 

2  C(K.J)=C(K,J)*TEMP 
IF(J.Eq.N)G0T01 
J=J+1 

G0T02 

1  R(K)=R(K)+TEMP 

J=1 

8  IF(K.EQ. J)G0T03 

TEMP=C(J,K) 

L=K 

5  C(J,L)=C(J,L)-C(K,L)*TEMP 
IF(L.Eq.N)G0T04 

L=L+1 

G0T05 

4  R(J)=R(J)-R(K)*TEMP 

3  IF(J.Eq.N)G0T06 
J=J+1 

G0T08 

6  IF(K.Eq.N)G0T07 
K=K+1 

G0T09 

7  RETURN 
END 
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* 

*^i:^:^^^^i^:^^i^^^t*^^*^i*****^:4i>|^:^^i^t^*********************************************** 

* 

♦This  subroutines  determines  if  the  angle  (phi  or  psi)  is  greater 
♦then  180  degrees.  This  is  required  because  of  how  the  arc 
♦cosine  function  is  computed. 

* 

subrout ine  t ruang (xc , yc , yr2 , xr2 , yr 1 , xr 1 , xdelt a , ydel t a , phi ) 
real  mref 

if (xr2 . eq . xrl) then 

if (xc . gt . xdelta)phi=360 . -phi 


else 


mref=(yr2-yrl)/(xr2-xrl) 
bref =yc-mref ♦xc 
ck=mref ♦xdelt a+bref 
if (mref . ge . 0 ) then 

if (ydelta.gt . ck)then 
phi=360 . -phi 
endif 


else 

if (ydelta . It . ck)then 
phi=360 .-phi 
endif 

endif 


endif 

end 
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APPENDIX  B: 

INSTRUCTIONS  FOR  USING  THE  YAW-CARD  REDUCTION  PROGRAM 
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Intentionally  left  blank. 
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The  following  instructions  and  example  are  provided  to  aid  in  the  use  of  the  i)iogram 
listed  in  Appendix  A. 

STEP  1.  Create  “RDTYPE.DAT”  file 

This  file  contains  information  about  the  projectile. 

The  file  has  the  following  lines: 

DESCRIB:  A  one-line  description  of  the  round 
FINSP:  The  diameter  (or  span)  of  the  fin  blades 
PLEN:  length  of  the  cylindrical  portion  of  the  projectile 
PRAD:  radius  of  the  projectile 
FINNO;  number  of  fin  blades 

Figure  B1  shows  an  illustration  of  the  above  variables. 

The  following  is  an  example  listing  of  a  “RDTYPE.DAT”  file; 

M900  projectile 
3. 

23.35 

.46 

6. 


PLEN 


1/2  FL- 


i _ i 

, _  ^ 

1 

'_J 

L 

— P' 

FINSP 

I - 

1 

: - ^ 

— V 

UflH 

1 

r 

FINNO  =  6 


Figure  Bl.  Projectile  Information  Required  to  Execute  Piogiaiii. 
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Note  that  the  variables  FINSP,  PLEN  and  PRAD,  must  hav'e  the  same  unit  of  measure. 
In  this  example,  the  unit  is  inches.  Also  note  that  FINNO  is  the  number  of  fin  blades,  not 
the  number  of  fin  blade  imprints  on  a  particular  yaw  card. 

A  new  “RDTYPE.DAT”  file  must  be  created  whenever  a  different  type  of  projectile  is 

use. 

STEP  2.  Digitizing  Pad  Layout 

The  layout  for  the  digitizing  pad  is  shown  in  Figure  B2. 


The  Areas  SI  and  S2  are  “switches”  required  by  the  program.  Si  signals  the  program 
that  all  the  usable  fin  blade  imprints  are  entered.  This  switch  is  required  because  there 
are  cases  when  the  round  (because  of  the  yawing  motion)  removes  the  material  on  the  yaw 
card  where  a  fin-blade  imprint  may  occur.  This  creates  a  missing  fin-blade  imprint.  The  S2 
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imprint  signals  the  program  that  all  the  data  from  one  yaw  card  are  entered.  This  switch  is 
provided  for  user  convenience.  It  allows  the  user  to  review  a  file  with  digitized  data  and  see 
where  the  data  from  one  yaw  card  ends  and  the  data  of  another  card  begins. 


STEP  3.  “Digitizing”  a  Yaw  Card 

The  points  that  need  to  be  digitized  are  shown  in  Figure  B3. 

The  digitized  data  goes  in  a  file  called  “DIGDAT.DAT”  in  the  following  order: 

Q1-Q2 

Reference  line  data.  Qly  >Q2y  and  Ql  must  be  entered  first. 

Q3-Q8 

Extreme  edges  of  the  fin  blade  imprints.  They  can  be  entered  in  any  order. 

Q9 

This  creates  a  flag  in  the  data  file  that  signals  the  program  that  all  the 
fin  imprints  are  entered. 


QIO 

Approximate  center  of  fin  imprints.  Required  to  start  iteration  process  to 
find  true  center  of  the  imprints. 

Qll 

Extreme  edge  of  the  yaw  imprint. 

Q12 

This  creates  a  flag  in  the  data  file  to  mark  the  end  of  the  data  for  this 
yaw  card. 
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The  following  is  an  example  of  a  “DIGDAT.DAT”  file: 


4.58343 

5.71357 

4.98731 

1.57662 

5.22843 

4.84069 

5.53013 

6.27565 

6.92620 

6.74218 

8.01478 

5.74371 

7.73103 

4.32088 

6.33700 

3.84737 

2.10548 

0.86783 

6.74561 

5.32114 

7.45783 

5.31157 

-1.34795 

-1.49891 

STEP  4.  Execution  of  Program 

With  the  files  “RDTYPE.DAT’  and  “DIGDAT.DAT”  created,  no  other  informa' ion 
is  required  by  the  program.  The  program  will  compute  8  ,Q  ,  /?,  0  and  i/>  when 
the  program  is  executed. 

A  sample  from  the  output  file  “OUTPUT.DAT”  is  shown  below; 

M900  projectile 

no.  delta  theta  alpha  beta  psi 

1  0.919  184.370  -0.070  -0.916  47.576 

where  no.  is  the  yaw  card  number. 

Although  only  one  yaw  card  result  is  shown  in  this  example,  several  yaw  cards’  digitized 
data  can  be  in  the  same  “DIGDAT.DAT”  file.  They  must  be  placed  in  that  file  in  the  proper 
sequence  (as  describe  previously  with  the  “Q”  points)  and  in  the  proper  yaw-card  order,  the 
first  yaw-card’s  data  in  the  beginning  of  the  file  and  the  last  yaw-card’s  data  at  the  end  of 
the  file. 
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Figure  B3.  Required  Digitized  Data  to  Execute  Program. 


Intentionally  left  blank. 
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APPENDIX  C: 

HOW  THE  £p  LENGTH  IS  DETERMINED 
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The  length  i-p  is  required  to  determine  8.  As  seen  in  Figure  Cl. 

4  =  (Cl) 

where  is  the  largest  distance  between  the  beginning  of  the  cylindrical  portion  of  tlie  rod 
and  the  location  on  the  projectile  where  the  fin  blades  are  at  their  largest  span.  .4  study 
was  conducted  to  determined  how  the  Q  length  is  best  determined.  It  was  assumed  that  the 
center  of  the  fin-blade  imprints  was  known. 

The  first  case  examined  is  a  projectile  that  has  two  fin  blades  in  the  yaw  plane.  This  is 
showm  in  Figure  Cl.  When  this  occurs,  the  length  of  the  cut  made  by  the  fin  blades  is  2L. 
The  location  of  the  vertical  center  of  the  fin  blade  imprints  is  L  and  as  stated  previously,  the 
L  value  is  assumed  known.  The  intersection  of  lines  //<-  and  Ipc  (the  centerline  of  the  flight 
body)  is  required  to  determine  Q. 


When  the  terms  are  combined, 


/  =  — ^ — sinS  —  (E  +  -D)smStanS 


(C12) 
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and 


M  =  l{A  +  B)-iE-\-l-D)tan6. 


(C13) 


Finally,  from  Equations  C2,  C3  and  C13 

1 


Q=-(A-^B). 


(C14) 


The  second  case  examined  is  a  projectile  with  a  two-bladed  fm,  where  the  blades  are  in 
the  plane  perpendicular  to  the  yaw  plane.  This  is  illustrated  in  Figure  C2.  In  this  case,  the 
center  of  the  fin  blade  imprints  is 


L  =  -{A  +  B)sinS. 


(C15) 


Again,  the  location  of  the  intersection  of  lines  I/c  and  Ipc  is  required  to  determine  Q.  The 
length  Q  is 


Q  = 


sinS ' 


therefore. 


Q=-(A  +  B}. 


(C16) 


(C17) 


The  results  of  these  two  examples  are  the  same.  It  is  assumed  that  if  the  center  of  the 
fin  imprints  {L  in  Figures  Cl  and  C2)  can  be  determined,  the  correct  Q  length  for  all  other 
cases  is 

(C18) 


Q  =  -(A  +  B). 
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LIST  OF  SYMBOLS  AND  ABBREVIATIONS 


Ai  and  Bi  are  two  of  the  direction  cosines  of  the  line  segment  that  has  Points  PI  and 

P2  (as  seen  in  Figure  3)  as  end  points.  and  B2  are  similarly  defined  for  the 
reference  line. 

6re/  y-intercept  of  the  line  that  has  a  slope  of  Are/  and  goes  through  the  center  of  the 
fin  blade  imprints. 

Cien  distance  between  the  extreme  edge  of  the  yaw-card  hole  and  the  center  of  the  fin 
blade  imprints 

(.p  projectile  length  required  to  determine  8.  See  Figure  Bl. 

it  largest  distance  between  the  beginning  of  the  cylindrical  portion  of  the  rod  and  the 

location  on  the  projectile  where  the  fins  are  at  their  largest  span.  See  Figure  Bl. 

Prii  the  “x”  value  of  the  P  point,  where  the  P  points  are  specific  points  on  a  yaw 

card.  See  Figure  3. 

Puy  the  “y”  value  of  the  P  point,  where  the  P  points  are  specific  points  on  a  yaw 

card.  See  Figure  3. 

Qn  points  on  the  yaw  card  that  must  be  digitized  in  order  to  execute  the  program  in 
Appendix  A. 

Q  it  —  ip'i  it  is  5FL  on  Figure  Bl. 

r  the  largest  radius  of  the  flight  body 

Rfin  radius  of  the  fin-span  circle.  It  is  FINSP  in  Figure  Bl. 


31 


Greek  Symbols 

a  pitch 

j3  yaw 

6  total  angle  of  attack 

Are/  slope  of  the  reference  line  in  the  digitizer’s  coordinate  system 

(f)  angle  between  the  reference  line  and  the  line  that  contains  Points  PI  and  P2  (see 

Figure  3).  The  angle  begins  at  the  reference  line  and  goes  in  a  clockwise  direction. 

(f)'  principle  value  of  the  arc  cosine  function 

4^  angle  between  the  /3-axis  and  the  first  fin  blade  imprint  clockwise  of  the  /3-axis 

0  angle  between  the  /3-axis  and  the  line  that  contains  Points  PI  and  P2  (see 

Figure  3).  The  angle  begins  at  the  /?-axis  and  goes  in  a  clockwise  direction. 

Oc  windsheild  angle 
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